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Abstract 



c I ■ The total cross section for top quark pair production close to threshold 

in e + e~ annihilation is investigated. Details are given about the calculation 

Q_i! at next-to-next-to-leading logarithmic order. The summation of logarithms 

leads to a convergent expansion for the normalization of the cross section, 
and small residual dependence on the subtraction parameter v. A detailed 
analysis of the residual v dependence is carried out. A conservative estimate 
for the remaining uncertainty in the normalization of the total cross section 
from QCD effects is Sa t t/a t t < ±3%. This makes precise extractions of the 
strong coupling and top width feasible, and further studies of electroweak 
effects mandatory. 



Typeset using REVTgX 



I. INTRODUCTION 

Detailed studies of the top quark are among the major projects of a future lepton pair 
collider. Of particular interest is the production of top quark pairs close to threshold, where 
the top quark velocity v is small and the usual perturbative expansion in terms of the strong 
coupling breaks down due to Coulomb singularities. The large top quark width T t ~ 1.5 GeV 
prohibits the production of toponium states and leads to a cross section that rises smoothly 
when top pair production becomes kinematically allowed. Furthermore, the large decay rate, 
T t ^> Aq CD , serves as an infrared cutoff, allowing the use of perturbative methods for the 
description of the non-relativistic top-antitop dynamics to a high degree of precision JTJ, Q] . 

In the past numerous studies have been carried out to assess the feasibility and precision 
for extracting various top quark properties from a threshold run |5|-[J. Recently, next-to- 
next-to-leading order (NNLO) QCD corrections to the total cross section were calculated 
using the concept of effective field theories |TD|-|TIJ . Surprisingly, the corrections were found 
to be as large as the next-to-leading order (NLO) QCD corrections. From the residual scale 
dependence in the NNLO result, the normalization of the cross section was estimated to have 
~ 20% uncertainty |T7] . It was concluded that top quark short-distance mass parameters can 



be determined with a precision of 200 MeV or better from the shape of the cross section, 
if so-called threshold mass parameters are employed fT7|, ITS]. However, the large NNLO 
QCD corrections to the normalization of the cross section seemed to jeopardize competitive 



measurements of the top width, strong top coupling [I8|, or the top Yukawa coupling in the 
case of a light Higgs. Moreover, the large NNLO corrections seemed to indicate that, despite 
the perturbative nature of the tt system, high precision computations are not feasible. 

A common feature of all the NNLO QCD calculations is that the running from the 
hard scale m t down to the non-relativistic scales which govern the dynamics of the top- 
antitop system was not taken into account. In other words, at NNLO potentially large 
QCD logarithms of ratios of the hard scale and the non-relativistic scales m t v and m t v 2 
were treated perturbatively. From weak decays of K and B mesons, it is known that a 
consistent summation of logarithms between mw and the low-energy hadronic scales can 
significantly change the magnitude of Wilson coefficients and the corresponding physical 
predictions (see for e.g. Ref. |1|). In the case of the top quarks near threshold all the scales 
are well separated, m t ^> m t v ^> m t v 2 ^> Aqcd- Since ln[m t /(m t v 2 )] ~ In (40) = 3.7, these 
logarithmic contributions can be sizeable (see for e.g. Ref. pT 



In this work details are presented for the computation in Ref. [21j where the impact 
of the summation of QCD logarithms of v on the photon induced total cross section were 
examined. We also add the vector and axial-vector contributions from Z exchange. In this 
framework the expansion for the normalized cross section R takes the form 

R= -^- = ^Ef-) fc E(« S ln ^ x { 1 (LL);« s ^(NLL);a2 5a ^ ; ^ (NNLL) j 5 (1) 
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where the indicated terms are of leading logarithmic (LL), next-to-leading logarithmic 
(NLL), and next-to- next-to-leading logarithmic (NNLL) order. The summation of loga- 
rithms can be performed using renormalization group equations in the framework of "veloc- 
ity NRQCD" (vNRQCD) |f22}-|2~5| , an effective theory for heavy quark pairs. The logarithmic 



corrections treated in this work are at NNLL order and fully include the known NNLO QCD 



corrections from Refs. [|lT)|-|T6|| . A complete summation requires knowledge of the Wilson co- 



efficients for potentials and non-relativistic production currents at this order. For the NNLL 
potentials the required Wilson coefficients have been calculated in Refs. |23| , |25| ,p6|l. The LL 
running of subleading ti production currents is also required and is given below. For the 
leading production current the Wilson coefficient is known at NLL 



g and NNLO f2^9 

For this current only partial NNLL results are known and an estimate is made for the un- 
certainty this induces in our result. The combination of known NNLL results significantly 
reduces the size and scale dependence of previous NNLO QCD predictions |L7 . 



We emphasize that we do not treat electroweak effects consistently at the same level 
as the QCD corrections, as implied by the counting T t ~ m t v 2 used in this paper. This 
would require the inclusion of non-factorizable corrections as well as the single- and non- 
resonant background contributions, which is beyond the scope of this work. Here we only 
include the top quark width through an imaginary mass term, which can implemented by 
the replacement 0, E = \fs — 2mt — > E + iT t into a calculation for stable quarks. For the 
total cross section this represents a consistent next-to-leading order treatment of electroweak 
effects |30|]. We stress, however, that with the small uncertainties in R resulting from our 
NNLL QCD results, a systematic treatment of higher order electroweak related effects is 
desirable. 

The outline of the paper is as follows. In section [TT] we review an example from inclusive 
non-leptonic i?-decays to emphasize some important differences between fixed order and 
renormalization group improved calculations. In section |I| we discuss the matching and 
running of Wilson coefficients for top quark production near threshold. In section [IV] we 



evaluate the time ordered products of currents at the endpoint of the renormalization group 
evolution. Pole mass results are given in section [Vj, while our implementation of the IS 
mass is presented in section [VI]. In section |V11| the NNLL IS mass results are presented 
and the remaining theoretical uncertainty is analyzed. The consequences of our results for 
measurements of the top width, a s (mz), and the Higgs top Yukawa coupling are discussed 



in section VIII, followed by conclusions in section IX. 



II. SUMMING LOGARITHMS 

To outline the steps in our renormalization group improved calculation it is instructive to 
draw a parallel with the well-known computation for inclusive non-leptonic B-meson decays 
via the quark decay b — > cud. This process involves several scales — mw, rrib, rn c and Aqcd, 
but only mw and rrib will be discussed here. The relevant weak decay Hamiltonian at the 
scale /i is 
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where Pl = (1 — 7s)/2. Hw is /^-independent and the \i dependence of Cj(/x) is compensated 
by the \i dependence of Oj(/i). The coefficients Ci^irnw) are determined by integrating out 
W bosons in the b — > cud amplitude at ji = mw 



Ci (/i = m w ) = 1 + O (a s (m w )) , C 2 (fi = m w ) = + O (a s (m w )) ■ 



(3) 



The coefficients can be computed using perturbation theory in a s (mw), which is valid as long 
as a s {mw) is small. There are no large logarithms in the matching calculation since we have 
chosen \x = mw- The decay amplitude can then be computed by taking the matrix elements 
of 0\ (/i = mw) in the 5-meson. The problem is that these matrix elements contain large 
logarithms of mw/m, b . Using the renormalization group these large logarithms can instead 
be moved into the coefficients C, and summed by scaling Oi and Ci from \i = mw to /x = m&. 
Note that in a fixed order expansion without the renormalization group large logarithms exist 
even for \x = m b . The renormalization group equation (RGE) for the coefficients is 



and is diagonalized by taking C± 
over the integration range in Eq. ( 
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G\ ± C%. Provided the strong coupling a s {p) is small 
, higher order terms in 7 and (3 can be neglected and 
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Using a s (mw) — 0.12 and a s (mb) = 0.22 the coefficients C± at /i = m& are 

C+(m b ) = 0.85, C-(m 6 ) = 1.37, 



(5) 



(6) 



B i Jd 4 xTH w (x) H w (0) 
product expansion gives plf 



a change of -15% and 37% from their values at mw- Finally, the inclusive B meson decay 
rate is computed from the imaginary part of the matrix element of a time-ordered product, 
w ( x ) Hw(0) B)]. At leading order in l/m b at the scale \x = m fe , the operator 
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where f{m 2 c /ml) is a phase-space factor. The QCD corrections to this resummed result are 
of order a s (m b ) with no large logarithm. 

In the above example the computation of the decay rate was divided into three parts: 

1. Compute the coefficients Cj (// = m w ) in a perturbation series in a s (m w ), which is 
valid as long as a s (mw) is small. 

2. Scale Cj (/i) from \x = mw to // = m b using the renormalization group, which can be 
done provided a s (/i) is small in the region m b < fj, < mw- This scaling sums terms 
of the form [a s hi{mi,/mw)] k ■ The terms neglected in the anomalous dimension are 
smaller than those retained by a s (fi) at all points in the integration region. 

3. Compute the decay rate using Hw renormalized at /1 = m b in a perturbation series in 
o> s {m b ) (and l/m b ), which is valid as long as a s {m b ) is small. 

Since a s (n) is monotonically increasing as \x decreases, all three steps above are valid as long 
as a s (m b ) is small. However, it is crucial to note that because of the renormalization group 
improvement we do not need to assume that 



a s (m w ) In 
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is a small expansion parameter. In contrast, a fixed order computation of the decay rate 
in powers of (x s (ji) contains logarithms of mw I H and m&//^, and is valid only if a s (/x) times 
these logarithms is small. In particular, one needs a s (m,b) ln(mfp / 'm 2 ) to be small. As a 
result, fixed order calculations have limited validity when processes involve widely separated 
scales. 

For non-relativistic top quarks the running spans an even larger range of energies than the 
above example, fi = m t down to fi = m t v 2 ~ 4GeV. Therefore, at the order that anomalous 
dimensions depend directly on the m t v 2 energy scale (NNLO), it is not too surprising that 
renormalization group improved perturbation theory is necessary to compute the it cross- 
section. Since E ~ 4 GeV, fixed order perturbation theory would be valid only if 
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m t 



E 2 



0.8 (9) 



were small. In contrast to Eq. (|9]), renormalization group improved perturbation theory only 
requires a s {m t v 2 ) to be small. Relative to 5-decays the radiative corrections for top quark 
production could be even larger, since factors of a s for the bound state do not always come 
with a 1/tt. This is true of both the constant and logarithmically enhanced terms. 

For the highest order computed, fixed order perturbation theory cannot distinguish the 
scale at which to evaluate a s . Mistaking an a s (m t ) for an a s (m t v 2 ) is a difference of a 
factor of two. Renormalization group improved perturbation theory distinguishes the a s 's 
that appear at each order. We note that in our NNLL order results, the summation of 
logarithms which involve the m t v 2 scale do tend to numerically dominate. 

III. EFFECTIVE THEORY 



For our calculation we will employ vNRQCD [p^ -p5|, an effective field theory which 
describes the non-relativistic dynamics of heavy quark pairs, where the non-relativistic scales 
mv (momentum) and mv 2 (energy) are larger than the hadronization scale Aqcd, m being 
the heavy quark mass. This effective theory contains a consistent power counting in v, 
so that a given Feynman diagram contributes to a single order in the v expansion. It 
also accounts for the fact that the momentum (soft) and the energy (ultrasoft) scales are 
correlated through the heavy quark equations of motions. Finally, vNRQCD allows for the 
summation of logarithms of ratios of the hard scale m and the non-relativistic scales using 
a velocity renormalization group. The most important momentum regions for heavy quark 
processes are those with (energy, momentum) that are hard (m, m), soft (mv, mv), potential 
(mv 2 ,mv), and ultrasoft (mv 2 ,mv 2 ) JEJ. The effective theory is formulated by including 
only those quark and gluonic degrees of freedom which can become on-shell for scales below 
m. The on-shell degrees of freedom are gluons and massless quarks with soft and ultrasoft 
energies and momenta, and heavy quarks with potential energies and momenta. All off-shell 
effects such as those from hard quarks and gluons, potential gluons, and soft quarks are 
accounted for by on-shell matching of vNRQCD to full QCD at the hard scale. 
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FIG. 1. Examples of potential (a), soft (b), and ultrasoft (c) interactions in vNRQCD. 



Only ultrasoft energies and momenta are treated as continuous variables, and ultrasoft 
momenta are separated from larger momenta by employing the multipole expansion. Soft 
energies and momenta appear as discrete indices for potential quarks and soft gluons and 
must be summed over, similar to the velocity index v^ for the static quark in HQET J31 . 
For instance, the first few terms that are bilinear in the quark field tp p are 

C(x) = X>J(*){tf° " ^ + £?}M*), (10) 

where m is the pole mass, and there is a similar equation for the anti-quark field Xv Sums 
over soft and potential energies and momenta in combination with ultrasoft integrations in 
loop integrals can be rewritten as continuous loop integrals of soft and potential energies and 
momenta. The vNRQCD Lagrangian also contains operators that describe potential-type 
4-quark interactions (Fig. 0a) originating from potential gluons and other off-shell modes, 
interactions of quarks with soft gluons (Fig. |I]b), and ultrasoft gluons (Fig. [TJc) , as well as 
self-interactions of soft and ultrasoft gluons. The potential and soft interactions depend 
non-locally on the soft indices, but are local in ultrasoft momenta. Gauge invariance in soft 
energies and momenta is recovered through reparameterization invariance. Loop integrals 
are regularized with dimensional regularization in d = 4 — 2e dimensions using the MS 
subtraction scheme. 

The crucial feature of vNRQCD is that the dimensionful MS subtraction scale for loop in- 
tegrations over soft and potential momenta is fis, whereas for loop integrations with ultrasoft 
momenta it is \Lxj. Within dimensional regularization the two different subtraction scales 



are necessary to properly recover the full d- dimensional non-relativistic phase space [22 



Powers of (f^uY and (fisY are determined uniquely by the interaction vertices |24j as shown 



in Fig. [I|. Furthermore, the two scales are correlated, fig = vnv and [ijj = mis 2 , where the 
dimensionless velocity subtraction point v takes the role of the MS scaling parameter. The 
correlation between \xs and \x\j is an inherent property of vNRQCD, since the energy and 
momentum scales are related through the quark equations of motion. For NRQED this cor- 
relation is necessary to reproduce through running the (lna) fc , k > 2 contributions in Lamb 
shifts, hyperfine splittings and the corrections to the ortho- and para-positronium decay 
widths [^3| , |34| . The renormalization group equations for Wilson coefficients are determined 
in the canonical way with the exception that they are formulated in terms of the subtraction 
velocity v. 

Lowering v to a scale Vq of order the quark velocity sums all logarithms of the soft and 
ultrasoft scales into the Wilson coefficients of the vNRQCD operators. At NNLL order this 
includes logarithms originating from ultrasoft radiation effects. The ti state is predominantly 

6 



Coulombic, so vq ~ CVa s — 0.15-0.20. Once v is lowered to vo, power counting shows that 
matrix elements with ultrasoft gluons no longer have to be taken into account at NNLL for 
the description of a heavy quark pair. Thus, at v = Vq the NNLL Green functions of a color 
singlet quark-antiquark system are described by the two-body Schrodinger equation with 
relativistic corrections. 

In the remainder of this section we introduce the Wilson coefficients which are matched 
at the hard scale iy = 1), and run from v — 1 to Vq. 



A. Potentials 



The potentials in the NNLL Schrodinger equation arise from the four-quark matrix el- 



ements of the potential-type operators |35| and from time-ordered products of operators 
describing interactions with soft gluons. The latter account, for example, for the radiative 
corrections to the Coulomb potential, and include the dependence of the Schrodinger po- 
tentials on logarithms of soft momenta. Adding these contributions, the effective order 1/v 
Coulomb potential is 
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where k = p—q is the momentum transfer. The first term in Eq. (ffl|) has a coefficient V^(^) 
that contains the summation of NNLL logarithms coming from soft and ultrasoft gluons and 
was determined in Ref. p6 . The second term contains the one- and two- loop corrections to 
the Coulomb potential |36|. In vNRQCD they arise from time-ordered products of the lowest 
order operators describing the interaction of quarks with soft gluons |26[ . The couplings in 
these interactions are simply a s (fis) = ot s {mv) and evolve with the QCD /5-function. 

At NNLL order a virtual photon or Z boson produces top quarks in 3 5*i and 3 Pi states 
through the vector and axial-vector currents respectively. We can simplify the spin depen- 
dence of the subleading 1/m 2 potentials by projecting onto the spin triplet channel. To do 
this we choose the scheme where traces over a 1 matrices are done in three dimensions (using 
S 2 = 2 for triplet states). 1 The order v° and v 1 potentials are0 
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V, 
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1 The difference between using three and d dimensional a matrices is simply a change in renor- 



malization scheme. A similar scheme dependence arises in chiral perturbation theory [37]. 



2 At NNLO the potentials V$, V r and Vk are equivalent to what was referred to as "Breit-Fermi" 
and "non-Abelian" potentials in some publications on the NNLO corrections (see e.g. |10-12|1). 



In vNRQCD the potential coefficients at v — 1 are obtained with on-shell matching at 
the hard scale. The potentials V$ and V r are of order a s and lead to terms in the cross section 
that are v 2 -suppressed. Their evolution therefore only needs to be taken into account in LL 



order and was determined in Ref. [23]. The potential V& arises only at order a s and also 



leads to terms in the cross section that are v -suppressed. Its evolution needs to be known 



at NLL and was determined in Ref. [5S| . At NNLL order the complete sum of potential-like 
interactions is 

V(p, q) = V c (p, q) + V 5 {p, q) + V r (p, q) + V k (p, q) . (13) 



Explicit formulae for the Wilson coefficients and constants in Eqs. flllD and (12) are collected 
in Appendix A. 

B. Currents 

To describe quark-antiquark production close to threshold at NNLL order we need to 
know the Wilson coefficient at NNLL order of the dimension three 3 S"i production cur- 
rent, and at LL order the dimension four 3 Pi and dimension five 3 Si currents. The vector 
production current is J p = c\O pA + c 2 Pt2 , where 

Pjl = ^cr(ia 2 )x- p , (14) 

O p2 = — i>p P 2 (r(ia 2 ) x*-p , 

and the relevant axial- vector current is J p = c 3 Pj3 , where 

O p ,3 = — V [tT,a--p]{ia 2 )X-p- (15) 

In Eqs. (|I^) and (^) p is the soft-momentum index of the quarks. The corresponding 
annihilation currents Pi i_3^ are obtained by complex conjugation. In the basis of operators 
we are using there is an additional dimension five vector current, 

1 . tt L,. ^ -P 2 



O pA = — 2 V [p{<r ■ p) - a- ij j (ia 2 ) x-p ■ (16) 

However it produces a Z)-wave quark-antiquark pair and therefore does not contribute at 
NNLL order. 

The matching condition at the hard scale for the Wilson coefficient c\[y = 1) is needed to 
order a 2 s . Its value is determined from matching the two-loop result for the quark-antiquark 
production amplitude close to threshold in full QCD to the corresponding amplitude in vN- 
RQCD and is given in Appendix A. Since the two-loop result for Ci(l) is scheme dependent, 
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FIG. 2. Graphs with an ultrasoft gluon with p ■ A couplings which contribute to the running 
of c 2 (^). 



our result is different from that obtained with the threshold expansion J3^j . The LL anoma- 
lous dimension for C\ is zero. The evolution of c\ for v < 1 at NLL order was determined 
analytically in Ref . |25| . Due to the length of the resulting expression this analytic result 
for v < 1 is not repeated here. At NNLL order the evolution of c\{v) requires the three loop 
anomalous dimension which is not completely known. In Appendix A we discuss the form 
of contributions to this anomalous dimension, and investigate their effect on Ci(u). Based 
on the contributions that are currently known, and a parametric estimate for the size of the 
unknown contributions, we find that the uncertainty in C\{y) from neglecting this running 
is likely to be at a level < 1%. 

The coefficients c 2 and C3 are needed at LL order. This involves tree level matching, 
together with one-loop running. At tree level expanding the full theory currents, tyt and 
t77 5 t gives 



c 2 (v = 1) = -1/6 . 



c$\y 



1. 



(17) 



The LL evolution of c 3 would be determined by divergent order a s v graphs, but no such 
diagrams exist so c^{v) = 1. For c 2 (z/) we need the divergent order a s v 2 graphs which are 
shown in Fig. 0. Together they give the anomalous dimension 
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Integrating Eq. ( |T8"D with the boundary condition in Eq. flT7| ) and the LL value of C\(v) gives 
the LL solution for c 2 (^)^ 
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A check of this result is obtained by expanding in a s (m) 



C2[y) = —z + -=— a s {m) In ( 

07T m 
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where the ln(z/ 2 ) term agrees with the ln(/x/m) term derived in Ref. p8[. (Note that the 



operator O p2 in Ref. |38[ was defined with the opposite sign. 



Eq. ( |19| ) was used in Ref. |Ell, but corrects a typo in Eq. (13) of this reference. 



C. Width and Nonperturbative Effects 

To incorporate the effect of the large top quark width we include in the effective La- 
grangian the operators 



V V 



£ = £vi^ p + E4^ tXp . (2i) 



In the Standard Model the dominant decay channel is t — > feW 74 " and gives a width of 
r t = 1.43 GeV (including one-loop electroweak |[j| and two-loop QCD J4(| effects). For 



power counting we will treat r t as a fixed number to be determined from experiment. The 
typical energy of a Coulombic top quark is E ~ mv 2 ~ 4 GeV, and we take Y t ~ E. Thus, 
the propagator for a single top (or antitop) with momentum (p°,p) is 

(22) 



p° - p 2 /(2m) + iTt/2 + it 



In our results the use of this propagator is equivalent to the replacement E — > E + iY t in 
ti Green functions from the stable quark calculation. For the total tt cross section in the 
Coulombic regime the operator in Eq. fl21"|) is sufficient to account for all electroweak effects 



at LO || and NLO |30|. The complete set of electroweak effects at NNLL order is currently 
unknown. 

The large top quark width plays a crucial role in suppressing the size of nonperturbative 
hadronic contributions governed by the scale Aqcd- Even at threshold (E = 0), the tt 
system has an effective energy set by the perturbative scale T t . For stable quarks, an 
operator product expansion fl4~i~l,[4"2l in powers of Aqqb/E can be used to incorporate non- 
perturbative contributions. Naively, the size of the first non-perturbative correction for tt 
production is therefore [Aqcb/{E + iT t )} 4 from the gluon condensate operator. However, 
because the momenta of the gluons, of order Aqcd, is much smaller than the momenta of the 
quarks mv ~ 25 GeV, their coupling to quarks is governed by a multipole expansion. The 
lowest order interaction is [ippP- A.ip p ], giving an additional v 2 suppression [fil^j . The size of 
the leading nonperturbative effect is therefore 5R np /R ~ Aq C - d v 2 /(E + iT t ) A i$ 10 -4 . This 
is consistent with the size of the gluon condensate corrections to R computed in Ref . [ 33 1 , 
and is more than two orders of magnitude smaller than the corrections considered here. 

IV. TOTAL CROSS SECTION 

In this section we compute the time ordered products of currents at a scale corresponding 
to the typical momentum and energy of the non-relativistic quarks, namely v ~ 0.15 — 0.20. 
In full QCD the expression for the total cross section a^t (e + e~ — > 7*, Z* — > qq) for quarks 
at center of mass energy y/s is 

^ot(s) = a pt [E^s)R v (s) + E a (s)R a (s)} , (23) 

where <7 p t = 47ra 2 /(3s). The vector and axial- vector i?-ratios are 
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where g = (y/s, 0) and j" (Jfy is the vector (axial- vector) current that produces a quark- 
antiquark pair. With both 7 and Z exchange the prefactors in Eq. ( f23"|) are 
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(26) 



Here Q/ is the charge for fermion /, T/ is the third component of weak isospin, and 6yy is 
the weak mixing angle. 

In vNRQCD at NNLL order the current correlators are replaced by the correlators of 
the non-relativistic currents P i, so that 
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Im 



Im 



cl(v)Ai(v,m, v) + 2c\{v)c 2 {v) A 2 {v,m, v) 
cl(is)A 3 (v,m,is) 



(27) 
(28) 



with 



A±(v, m, v) = i y^ / d 
p,p' 



4 x e^ x ( 



A 2 (v,m,v) = - J2]d A xe^ x (0 
p,p' 



A 3 (v,m, v 



J2 [d 4 x e iq - x ( 
p,p' 



TO p , 1 (x)O^ 1 (0) 
T [O pil (x)0;, )2 (0) + O p , 2 (x)O pv (0) 
TO p , 3 (x)O p , i3 (0) 0). 



(29) 



Here q = (y/s — 2m, 0), and the current operators O pi are defined in Eqs. ([14]) and (|T5|). 
The correlators Ai are functions of the quark mass m, velocity renormalization scale v, and 
the effective velocity 



s — 2m + iT t \ 2 



m 



(30) 



Here m is the pole mass. For our calculation evaluating matrix elements at the low-scale 
corresponds to computing the Ai's with v of the order a typical Coulombic top quark's 
velocity. The correlator A 2 can be related to A\ using the quark equation of motion giving 



A 2 (v, m, v) = v 2 Ai(v, m, v) 



(31) 
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Performing the spin traces in 3-dimensions the correlators A\ and .A3 are determined by the 
zero-distance coordinate space Green functions, 

Ai(v,m,u) = dN cf xf J D n pD n p' G v>mjV (p,p') , 

12 N r 

A 3 (v,m,is) = -^tfs J D n pD n p> '{p-p')G m (p lP ') . (32) 

Here n = d—1 = 3 — 2e, we use D n k = e elE (An)~~ e d n k/(2ir) n which converts /is (and v) 
from the MS to the MS scheme, and N c = 3 is the number of colors. The dependence of the 
A3 prefactor on n comes from taking the trace over a matrices in the P-wave current in 3 
dimensions, and then projecting the dot product p ■ p' back into n-dimensions. 

To determine the correlators A\ and A3 at NNLL order we use a combination of nu- 
merical and analytic calculations. Ultraviolet divergent contributions are computed with 
dimensional regularization in the MS scheme, which is required for consistency with the 
scheme used to compute the matching and running of the Wilson coefficients. Our strategy 



for the analytic calculations for these potentials follows that in Ref. [10|. In momentum space 
the NNLL Schrodinger equation for the Green function for a color singlet quark- ant iquark 
pair is 



G,,„,AP, p') + / D"k £• V(p, h) G,, m ,„(fc, p') = (21)" «<»>(p - p') , (33) 



m 4m d 

where E = mv 2 and V(p, k) is given in Eq. (|i~3"D. 

We start with the vector correlator A\. The lowest order S'-wave Coulomb Green function 
at zero distances reads 

Tfl ( r / — 2 V \ 1 / % d ' 

G°(a,v,m,u) — — < iv — a In I — - +hi2 + j E + ip [1 — — - 

47r L L V v ) 2 \ 2v/ 

m 2 a 1 ,„,. 

+ "J— T" 34 

where a = — V c (z/)/(47r). Eq. (Q) is a standard result |4^| , except for the scheme dependent 



constants entering with the divergence. This term was computed from the two-loop graph 
with a single V c insertion, and Eq. (J34j) agrees with the MS result obtained in Ref. p6| . The 
1/e pole is canceled by a counterterm for the time ordered product of currents. An identical 
divergence structure appears with the full V c potential in Eq. ([11]). Since this divergence and 
corresponding constants appear only in the real part of G°, for the purpose of computing 
Im[Ai] it does not matter whether or not they are subtracted in MS. Therefore, we determine 
the contribution from the Coulomb potential by the exact solution of the equation 

^ ~ E]G:^ u (p,p') + J -^V c (p,k)G c v ^ v (k,p') = (27r) 3 ^(p-p'), (35) 

using numerical techniques developed in Refs. ||£|]. We will denote by G c (a, v, m, v) the so- 
lution for the zero-distance coordinate space Green function obtained from G c (v , m, v){p,p')- 
On the other hand, the first order perturbative corrections to Ai(v,m, v) from the po- 
tentials \4, Vs, and V r , and the kinetic energy correction, p 4 /(4m 3 ), have ultraviolet sub- 
divergences which must be subtracted in MS. For the first order corrections to the zero- 
distance Green function in dimensional regularization from the kinetic energy correction and 
the potentials 14, V$, V r (leaving out the Wilson coefficients) we find: 
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4 V4e 2 



9 v d 



f-|l + a- + i -IG°(a,,,m,,) 



(36) 
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z a 

2^ 



(37) 
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l..l-^)-5+b.2+T* + *(l- ,,. 



SG r (a, v, m, v) 



256vr 2 e ' 
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(38) 



m 



167T 2 
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iv — a 



| ll( _)- 5+ ln2 +7 , + «(l-- 



+ 



771 
167 2 



V + 



0,- 



1 1 



4 V4e 2 



w 2 d 
47r 9a 



G (a, t>,m, v) . 



(39) 



In deriving these equations we have included the counterterms generated by renormalizing 
the O p i current |22 at NLL order. These counterterm graphs are sufficient to cancel all 
sub divergences. The remaining overall divergences in Eqs. ( |36| ) to ( |39"D are of the form 1/e 
and f 2 /e and are canceled by vacuum counterterms for the current correlator, and can be 
dropped. The final renormalized expression for the NNLL vector correlator A\ is then 



Ai(v, m, v) = 6 N c 



G c (a, v, m, v) + (v^ s) (u) + 2V s (s) (i/)) 5G s (a, v, m, v) 



+ VJ: s \v)6G r (a,v,m,v)+ V ( k s) (u) 5G k {a,v,m,v) + o"G kin (a,t;,m, v) 



(40) 



For the computation of the axial- vector correlator A3 at NNLL we only need to consider 
the leading order Green function because the axial-vector current is already suppressed 
by one power of the velocity. As A3 describes P-wave quark- ant iquark production it is 
proportional to the I = 1 component of the Coulomb Green function at zero distance, 



m 



G 1 (a,v,m,iy) =-- - — <; / r" a r" 

4:71 
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7 T / io 
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4 V 2v 



9 a 



(41) 



The divergences and constants were determined by explicitly computing the two-loop and 
four- loop potential graphs in the MS scheme (with one and three V c insertions respectively). 
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FIG. 3. Results for the vector current i?-ratio in Eq. ([27]) with a fixed pole mass, rrit = 175 GeV. 
In (a) we reproduce the NNLO pole mass results (with matching values for Wilson coefficients). In 
(b) we display our NNLL pole mass results (with running coefficients). In both plots the dotted, 
dot-dashed, dashed, and solid curves correspond to v = 0.1,0.15,0.2, and 0.4 respectively. 



As before, the result contains overall divergences of the form 1/e and v 2 /e that are canceled 
by vacuum counterterms. The final result for the axial-vector correlator A3 is 



Az(y,m,v) 



4N r . 



m z 



G 1 (a, v,m, u) . 



(42) 



The (fixed order) NLO and NNLO corrections to .A3 are known |L6], 47]. Since they contribute 



.7,2 



to o"7ot beyond NNLL order, they are not included in our analysis. 

We emphasize again that the use of Eq. (^) in determining R v,a does not provide a 
consistent treatment of the top quark width at NNLL order. This can be seen from the 
presence of //-dependent terms proportional to a s Y t /m t ln(z//a) in the imaginary part of 
Eq. ([H]) JTOf , and from a scheme dependent term proportional to a s Y t /m t in the imaginary 
part of Eq. (ffTF). These terms are energy independent and parametrically of NNLL order. 



Conceptually, they indicate that additional renormalization and matching is necessary to 
achieve a consistent treatment of all electroweak effects at this order. The issue of electroweak 
effects in the total cross section will be addressed elsewhere. 

V. POLE MASS RESULTS 



Let us briefly examine the results for R" and R a in the pole mass scheme. In Fig. §a 
the NNLO photon-induced cross section Q\ R v is displayed over the cm. energy yfs for 
a s (m z ) = 0.118, pole mass m t = 175 GeV, T t = 1.43 GeV, and for v = 0.1,0.15,0.2 
and 0.4. For the strong coupling four-loop running is employed and all light quark flavors 
(rif = 5) are taken massless. Following Ref. |l7j] we have defined NNLO by taking v = 1 
in the current Wilson coefficients Ci^v) and in the ln(— ivjv) terms in the Green functions 
SG km ' k,s ' r . A lower scale /z so f t = m t v is used for factors of a s that appear in both the Green 
functions and in Vi{y — 1). At NNLO the size of the corrections and the dependence on v 



are found to be similar to previous NNLO calculations ||17|| , indicating that our calculational 
strategy does not change the fixed order results. 
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In Fig. [3|b we show the renormalization group improved NNLL pole mass results for 
the same parameters and range of v used in Fig. |3|a. From the physical point of view the 
appropriate choice of the subtraction velocity parameter v is around Cpa a ~ 0.15-0.2. With 
this choice all large logarithms of ratios of the the hard and the non-relativistic scales are 
summed in the Wilson coefficients, and the matrix elements of vNRQCD operators are free 
of large logarithmic terms. To be conservative we show a larger range of values, v = 0.1-0.4, 
in Fig. |3|. In both the NNLO and NNLL results the dependence of the peak location on v is 
a reflection of the use of the pole mass and is considerably reduced if threshold masses are 
employed [|T7J as we do in the next section. From Fig. ^ we already see that the dependence 
on v of the normalization of the NNLL results is significantly reduced in comparison to the 
NNLO results. This is discussed in more detail in the following section. 

An issue we would like to discuss here concerns the flattening of the peak that is observed 
for the NNLL v = 0.1 curve in Fig. |3|b. This flattening turns out to be specific to the use of 
the pole mass together with our calculational strategy. To see this consider the contribution 
from the n = 1, 3 S'i bound state to R v . Up to trivial factors its contribution to the cross 
section based on Eq. (|40"D reads 



R v n= i oc Im 



«0)| s 



Er-E-iYi 



+ 



«o)l; 



SE m \y° 



E°-E-iT t (E°-E-iT t y 



(43) 



Here |\l/ c (0)| 2 and E c are the square of the wave function at the origin and the binding energy 
of the ground state obtained from Eq. (J3~5f), and 5|^(0)|^ and 5E m are the corresponding 
perturbative corrections from the potentials Vk, Vg, V r and the kinetic energy correction 
p 4 /4m 3 . The terms |\l/ (0)| 2 = m 3 a 3 /(87m 3 ) and E° = —m t a 2 /4 are the LL squared wave 
function at the origin and binding energy, respectively. Explicit formulas for E c and 5E m in- 
cluding the summation of logarithms of v have been presented in Ref . [26] . If the unexpanded 
NNLL Schrodinger Equation in Eq. 



is used, then Eq. ( fi3|) becomes 



K=i ex Im 



|# c (0)| 2 + c^ 



E c + 8E m -E-iT t 



(44) 



For E « E c and SE m > 30MeV the relative difference between Eq. (§4j) and fl43|) is of 
order -SE m /T t , because E° - E c « 1 GeV is large. For v = (0.1,0.15,0.2,0.3,0.4) and 
using the same parameters as in Fig ||a, we obtain 5E m = (143,27, —10, —33, —40) MeV. 
Thus for v = 0.1 our calculational method leads to negative contributions of about 10% at 
energies around the nominal peak position with respect to a calculational strategy which 
instead involves the contribution in Eq. (|44j). This is observed in Fig. |3|b for v = 0.1. We 
have checked that one can cure this flattening by summing the perturbative contributions 
in Eq. ( fl"3|) into the single energy denominator of Eq. (0) (see also |13], |16]] ) . However, this 
procedure sums terms that are subleading according to the power counting. 

Instead, we note that this effect would be small if E° — E c ^C 1 GeV. Thus, it can be 
regarded as a direct consequence of the pole mass scheme where this energy difference is 
large due to an infrared renormalon. In the subsequent section we will make use of the IS 
mass scheme 



1581 



together with our calculation strategy. In this case the flattening effect 
does not occur because properly introducing the IS mass makes the analogue of the third 
term in Eq. (031) explicitly zero. 
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FIG. 4. Pole mass results for the axial- vector current i?-ratio in Eq. (^7]). The dotted, dashed, 
dot-dashed, and solid curves correspond to v = 0.1,0.15,0.2, and 0.4 respectively. 



In Fig. |] the NNLL order axial-vector current induced cross section R a is displayed 
over the cm. energy y/s for m t = 175 GeV, a s (mz) = 0.118 and T t = 1.43 GeV for 
v = 0.1, 0.15, 0.2 and 0.4. As expected from the power counting, R a is at the percent level 
and suppressed by two orders of magnitude with respect to R v at cm. energies close to 
the IS peak and continuously grows for larger energies. No peak-like structure is visible 
since the P-wave wavefunctions are zero at the origin. For yfs < 348 GeV it is conspicuous 
that R a becomes slightly negative. This unphysical behavior of R a arises because we have 
not included a complete treatment of top width effects at NNLO. In particular there are 
short- distance corrections which in the non-relativistic theory correct for the true physical 
top phase space [T5| . Integration over the true top phase space, of course, leads to a positive 
cross section. As the difference between the approximate and true phase space integration 
involves only hard momenta, it corresponds to a short-distance correction. We note that 



unlike MS, in a hard cutoff scheme the cross section remains positive [H| because there is no 
way to generate a negative result by cutting off the integration region of a squared matrix 
element. In any case, like other known sub leading width effects these corrections are quite 
small, of order T t /mt ~ 0.01 for both R" and R a |15 . 



VI. THE IS TOP QUARK MASS 



Although the heavy quark pole mass is IR- finite [[|^,|50| and convenient to use in bound 
state calculations, it suffers from a linear sensitivity to small momenta, which leads to 
an ambiguity of order Aqcd associated to the notion of the pole mass |5Tfl . In practical 



bound state calculations based on the Schrodinger equation this ambiguity turns up as large 
r-independent perturbative corrections in the coordinate space Coulomb potential V c (r). 
These corrections grow asymptotically like -/is«"n! H52| . Because the top quark mass 
can be determined from the threshold scan at a future lepton collider with experimental 



uncertainties smaller than Aqcd [jig] , it is advantageous to abandon the pole mass definition 
and replace it with a so-called "threshold mass" . 

Threshold masses are short- distance mass definitions, i.e. their ambiguity is parametri- 
cally smaller than Aqcd? and the expansion relating them to the MS mass converges very 



16 



well. They remove the large perturbative corrections in a systematic way order by order 
in perturbation theory based on the fact that the total static energy .E^at — 2m t + V c (r) 
is free of an ambiguity of order Aqcd @0|. At the same time, threshold masses differ 
from the pole mass by corrections that are parametrically of order ma 2 and comply with the 
non-relativistic power-counting ]54| (see Eq. (|33"D, where all leading order terms are of order 
ma 2 ). In other words, threshold masses account for the fact that the quark is off-shell only 
by a small amount of the order v 2 ~ a 2 <C 1. In this work we use the IS top quark mass 



definition | I5| , which is defined as half of the mass of a perturbative contribution to the 
fictitious stable 3 Si toponium ground state. In recent NNLO calculations of the threshold 
cross section alternative threshold mass definitions have also been used such as the "kinetic 



mass 



55] , and the "potential subtracted mass" [54]] and variants of it |56|. (See Ref. [17 



for a comparison of threshold mass definitions.) 

The relation between the IS mass and the pole mass can be obtained from the 3 S\ ground 
state {n = 1) solution of Eq. (B3) and reads, 
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For convenience, we have displayed the LL, NLL and NNLL contributions separately. The 
term A^ NLL denotes the NNLL corrections from the Coulomb potential V c , and A^ NLL the 
ones from the potentials V$, V r and 14 and from the kinetic energy correction p 4 /4m 3 . For 
A NLL and A^ LL the difference between using a and Cf <y s (Ml s i>) only enters beyond NNLL. 
Using the non-relativistic power-counting, which is required to correctly eliminate the pole 
mass, the inverse of Eq. (051) reads, 



mt = M. 



is 



{l + [A L 



"A* 



+ 



(A LL )^ + A™ LL + A 



NNLL 

m 



(47) 



where the LL, NLL and NNLL corrections are each given in brackets. It is mandatory to 
use relation ( 477 ) at LL, NLL and NNLL for the LL, NLL and NNLL cross section at the 
same low-scale subtraction velocity v in order to guarantee the proper cancellation of the 
large corrections associated with the pole mass ambiguity. 
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In all of the Wilson coefficients for the potentials and currents the replacement mt — > M 4 1S 
suffices at NNLL order. This is because m t and M 4 1S differ by an order v 2 amount. Expanding 
the small logarithm, ln(l + All), gives 

a s (m t ) = a s (Ml s ) - ^a 2 (M t ls )A LL + . . . . (48) 

The expansions for a s (m t u) and a s {m t u 2 ) in terms of a s (M f ls z/) and a s (M 4 ls z/ 2 ) are analo- 
gous. Since A LL ~ v 2 the second term in Eq. ( f4"8"D gives a v 3 correction, and in the cross 
section enters at N 3 LL (or higher). 

Implementing the replacement in Eq. (]47| ) in the low energy Green functions is more 
involved. We implement the additional corrections in the Schrodinger equation in Eq. ( |33| ) 



following the calculational strategy used in Sec. |TV]. The corrections in Eq. ([47]) that arise 
from the Coulomb potential, A LL , A NLL and A^ NLL , are implemented through Eq. (pop, which 
we solve exactly. This accounts for the corrections arising from the definition of the energy 
E = ^/i — 2m t and the corrections coming from the leading p 2 /m t kinetic energy term. 
(For the latter only the LL correction A LL is needed.) The term A™ LL from the energy 
definition is treated perturbatively and leads to an additional correction to the correlator 
Ai that reads, 

8A1 (v, M t ls , v) = 6 N c 5G ls (a } v, M t is , v) , 

A NNLL J 

SG ls (a, v, M, 1S , v) = - =2— -f- G°(a, v, M, 1S , u) . (49) 

v dv 

5G 1S is an order v 2 correction since v~ l d/dv ~ v ~ 2 compensates A^ NLL ~ v 4 . Finally, in 
all f 2 -suppressed Green functions (Sec. [TV], Eqs. (|36|)-(|39|) and ^.2,3), the LL correction A LL 
needs to be included using 

= { w ) (50) 

rather than Eq. (|3(]) as the definition of the velocity, and the pole mass mt has to be replaced 
by M} everywhere else. This procedure correctly implements the IS mass scheme at NNLL 
order. 

In the IS mass scheme multiple bound state energy poles in A\ caused by the perturbative 
treatment of the potentials Vs, r ,k and the kinetic energy correction only exist for higher radial 
excitations (n > 2), but not for the ground state (n = 1) that is responsible for the IS peak 
visible in the lineshape. This can be seen from Eq. (f4*9"|), which exactly cancels the ground 
state double pole terms from Eqs. (J36l)-(|39l). The corrections that arise from removing 
the double pole terms for n > 2 by summing the A™ LL energy shifts into single energy 
denominators are found to be at the per-mille level and negligible. 

In non-threshold processes the MS top quark mass is typically a more convenient mass 
definition to use than a threshold mass. The IS mass is a scale-independent quantity, and 
its relation to the MS mass can be determined using the upsilon expansion [^SJ at /i = mt 



[y = 1). Explicit formulae for the determination of the top MS mass from the top IS mass 
can be found in Ref. 
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FIG. 5. Comparison of Q\R V with fixed M t is mass for the fixed order and resummed expansions. 
The dotted, dashed, and solid curves in a) are LO, NLO, and NNLO, and in b) are LL, NLL, and 
NNLL order. For each order four curves are plotted for v = 0.1, 0.125, 0.2, and 0.4. 
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TABLE I. Numerical values of Q^R" which appear in the NNLL results in Fig. |5p. 

VII. DISCUSSION 

In this section we carry out a detailed analysis of R v and R a in the IS mass scheme with 
the main emphasis on assessing the remaining theoretical uncertainties in our computation. 
In Fig. [5] we have displayed results for Q\R V over the cm. energy ^fs for M f is = 175 GeV, 
®-s{ m z) — 0.118 and T t = 1.43 GeV. For the strong coupling four-loop running is employed 
and all light quark flavors (rif = 5) are taken massless. Fig. |5|a shows results at LO (dotted 
blue lines), NLO (dashed green lines) and NNLO (solid red lines), while Fig. [5]b shows the 
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FIG. 6. Using the M t is mass a) shows the LL result for R a , and b) shows the total NNLL cross 
section in pb. The dotted, dashed, dot-dashed, and solid curves correspond to v = 0.1, 0.15, 0.2, 
and 0.4. 



resummed results at LL (dotted blue lines), NLL (dashed green lines) and NNLL (solid red 
lines) order. At each order four curves are displayed corresponding to v — 0.1, 0.125, 0.2, 
and 0.4. At sfs = 350 GeV the LO, LL, and NNLO upper through lower lines correspond 
monotonically to v — 0.1 through v = 0.4. In contrast, the upper through lower NLO lines 
correspond to v = 0.4 through v — 0.1. At NLL the lower and upper lines at ^/s = 350 GeV 
correspond to v = 0.1 and v = 0.2 (whereas v = 0.4 lies in between). Finally, at NNLL the 
lower line is v = 0.1, while the upper line includes both v = 0.2 and v = 0.4 at the displayed 
resolution. Our NNLO results agree well with previous NNLO analyses, see Ref. |L7| for a 
synopsis of previous NNLO results. In Table | we have summarized the NNLL values for 
Q\R V for yfs = 347, 350 and 353 GeV with the same parameters as used in the figure. The 
values v = 0.125 and v = 0.275 are included in the Table since at ^/s = 353 the value of 
QfR v is minimized near v = 0.125, while for ^ = 350 its value is maximized near v = 0.275. 

For R v the LL cross section is equivalent to the LO one because the one-loop anomalous 
dimension of c\ vanishes. The NLL cross section is obtained by eliminating from Eq. (f40|) 
the corrections from the 1/m-suppressed potentials, the kinetic energy term (p 4 /4m 3 ) and 
the operators V 2 and O p j,. For c\ the NLL evolution and the one- loop matching condition 
is employed. For the Coulomb potential, the coupling V c {y) is set equal to — 47rCVa s (M t ls z/), 
and only the one-loop corrections proportional to [a s (M 4 ls ^)] 2 have to be taken into account. 

In Fig. ||a we have displayed R a at LL for M t is = 175 GeV, a s (mz) = 0.118 and 
T t = 1.43 GeV and v = 0.1, 0.15, 0.2, 0.4. The lower through upper lines for yfs < 352 GeV 
correspond monotonically to v — 0.1 through v = 0.4. With the IS mass, R a dips below 
zero for a/s < 347 GeV. This dip has the same origin as for the pole mass and is discussed 
in section 0. Finally, in Fig. ^b we have displayed our total NNLL result for the 7 + Z 
induced cross section, using the same parameters as in Fig. |B|a. 

Comparing Figs. |5|a and [|b the variation of the normalization of R" with v at NNLL 
order is considerably reduced in comparison to the NNLO result. Equally important, the 
size of the NNLL corrections is considerably smaller than the size of the corresponding 
NNLO corrections. The latter indicates that at the low scale the perturbative corrections 
to the threshold tt cross section are converging. In a^ t the remaining variation of R" with 



20 



0.02 



0.00 






-0.02 



-0.04 



-\- 



0.1 



^ = 347 GeV 

I I I I I L 

0.2 0.3 

V 



0.4 






i i i | i i i | 


- 


0.02 


A . - 


- " ' _ 


ft; 


\ — - - - ' 


- 


O) o.oo 

< 

-0.02 
nnd 


fs = 353 GeV 

\ i i i i i i i i i 


/ ' 

/ 



0.1 



0.2 0.3 

V 



0.4 



FIG. 7. Relative scale dependence of various NNLL contributions to Q\R V . The contributions 
in Eq. (f40|) were divided into those from G c vmv (dashed lines) , 5G k (dotted lines) , and the sum of 
fiQ8,r,km,is (d t-dashed lines), while the solid lines denote the sum of these terms. 



v is still much larger than the variation of R a . Therefore, to estimate the remaining overall 
uncertainty in a^ t the variation of R a can be neglected, and for the rest of this section we 
discuss only R v . 

First we wish to explore the origin of the reduced v dependence. In Fig. [7|a,b,c Q\R V 
(solid red line) is displayed as a function of v for yfs = 347, 350, 353 GeV, using the choice of 
parameters of Fig. |5]. We have also displayed separately, the contributions to Q^R" coming 
from G c (dashed line), 5G k (dotted line) and the sum of SG 5,r ' kin ' ls (dot-dashed line). For all 
curves in Fig. |7| the respective values at v = 0.2 have been subtracted. The results show that 
for any cm. energy the ^-variation of the individual contributions to R v is much stronger 
than the ^-variation of the sum. As an example, for yfs = 350 and v > 0.15, one finds 
that the variation of the contributions from G c and 5G k cancels almost entirely, whereas 
the contribution from the sum SG s,r,kin ' ls is nearly ^-independent. On the other hand, for 
\/i = 350 and v < 0.2 the ^-variation is dominated by the contribution from 5G k which 
rapidly decreases, whereas the contributions from G c and SG s,r,km ' ls are small. The situation 
is quite similar for other energies, and results directly from the evolution equation of the 
Wilson coefficients. These evolution equations constitute a system of coupled differential 
equations, where the values of the Wilson coefficients become related to each other for 
v < 1. This is different from the fixed order NNLO calculations carried out in previous 
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FIG. 8. Position of the peak in the cross-section versus v at NNLO (dashed) and NNLL (solid). 
The vertical line at v = 0.15 is a physically motivated endpoint for the running. 



works [|I0| -|T6||, where the proper anomalous dimensions of the potentials and currents were 
not taken into account, and where the Wilson coefficients at the soft scale were obtained from 
simply evaluating their hard matching conditions at the soft scale. In particular, in the fixed 
order approach, there is no numerical cancellation between the scale-dependence that arises 
from the corrections to the Coulomb potential V c and the corrections from the potential Vk 
(see Eqs. (O) and (|12"D). In addition, in the fixed order approach, the contribution from 
Vk becomes quite large at low scales because it is proportional to the square of the strong 
coupling. In fact, in previous fixed order NNLO calculations Vk was identified as one of 
the major sources of the large positive corrections to the normalization of the total cross 
section. Using the VRG, on the other hand, the Wilson coefficient for Vk decreases by almost 
an order of magnitude between v — 1 and v ~ 0.2, and even vanishes for v » 0.17 [25|. 
This is one of the main reasons for the small size of the NNLL corrections. From the 
previous considerations we conclude that the cancellation of the ^-dependence and the size 
of the corrections that arises in the sum the various contributions to R v is a genuine and 
systematic QCD effect, and, in particular, that it is justified to take the overall //-variation 
of R v as a tool to estimate the remaining theoretical uncertainties. 

Let us now consider the range of v that should be considered for a suitable error estimate. 
As mentioned before, from the physical point of view the appropriate choice of the subtrac- 
tion velocity parameter v is around a s ~ 0.15-0.2. With this choice all large logarithms of 
ratios of the the hard scale m t and the non-relativistic scales mtv and m t v 2 are summed 
into the Wilson coefficients, whereas the matrix elements of vNRQCD operators are free of 
such logarithmic terms. From Fig. |6|b we see that if the lower limit for v is taken to be the 
endpoint of the running, v ~ 0.15, then there is really only a very small uncertainty from 
the residual scale dependence. To be conservative we feel that the v = 0.1 curve should be 
included. The curves in Fig. |7| show that our results become unstable for v < 0.1. This is 
a consequence of the dependence of the VRG equations on the ultrasoft scale nu = mtu 2 , 
which becomes of the order 1 GeV for v < 0.1. Thus, for v < 0.1 the perturbation expansion 
is apparently unreliable. However, this is just an artifact of choosing m t z/ 2 much smaller 
than the physical energy scales in it, so that a s {m t v 2 ) is large. The same feature can also 
be observed in Fig. §, where the position of the IS peak is displayed as a function of v for 
the parameter choice used previously. For v > 0.15 the peak position is stable, while for 
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v < 0.1 the peak position increases rapidly to higher energies. This is in contrast to the 
expectations in the IS mass scheme, where the peak position should be stable by construc- 
tion. The behavior of the peak position for v < 0.1 is a consequence of our perturbative 
treatment of the corrections that arise from the 1/m-suppressed potentials and the kinetic 
energy correction. For v < 0.1 these corrections become large and render the perturbative 
treatment unreliable. Thus, we conclude that the range v < 0.1 should not be used for the 
estimate of the remaining theoretical uncertainties. 

As far as values of v larger than 0.2 are concerned the curves in Figs. |5|-|8] show that 
our results are remarkably stable up to quite large values of the subtraction velocity, even 
beyond v = 0.4. For the estimate of the remaining theoretical uncertainty of the total cross 
section it is therefore not mandatory to fix any strict upper bound of v. In Table, [j] we have 
shown results for v = 0.1, 0.125, 0.2, 0.275, and 0.4. We consider the small value, v = 0.1, 
as a conservative lower bound, which is set by physical considerations, and below which no 
sensible perturbative treatment is possible. A less conservative approach would have been 
to take the lower bound to be set by the typical top quark velocity, v ~ 0.15 — 0.2. The 
value v = 0.275, on the other hand, has been included because for this choice the peak cross 
section reaches its maximum as a function of v. We take the size of the resulting v- variation 
of the cross section at the IS peak, 

— = ±3% (51) 

Ctt 

as the remaining uncertainty inherent to our renormalization group improved calculation. 
This is an order of magnitude smaller than the uncertainty that was associated to previous 
fixed order NNLO QCD calculations of the cross section. Note that if v = 0.125 was instead 
used as the lower bound for the running then the v- variation is Sau/crtt = ± 1.5%. In fact, 
the small ^-dependence and the small size of the NNLL order corrections observed in our 
calculation of the cross section comply with the expectations from power counting and the 
natural notion that the non-relativistic ti dynamics should be calculable to a high degree of 
precision [fllfl. 

It is instructive to compare the error assigned to the cross section based on the v- variation 
with the size of corrections from terms that are not yet included in our NNLL result. Firstly, 
at NNLL the complete running of C\[y) was not included because the three loop anomalous 
dimension for C\[y) is not yet completely known. (We recall that already at NLL order 
\c\[y ~ &s)} 2 contains the sizeable negative (fixed order) N 3 LO normalization correction 
proportional to o^ln 2 a s that was determined in Ref. [^Ofl , see Ref. [p5fl .) Contributions to 
the three loop anomalous dimension are discussed in Appendix |A2|, and based on dimen- 



sional analysis and the known contributions to this anomalous dimension we estimate that 
the uncertainty in neglecting it is ;$ 2%. It is also useful to consider the size of some correc- 
tions from beyond NNLL that can be determined easily or have been obtained earlier. We 
emphasize that all such corrections belong to matrix elements of the effective theory, i.e. it 
is crucial to employ the proper choice for the renormalization scale n$ or /i^/ depending on 
whether the corrections come from the soft (or potential) or ultrasoft momentum regime. 
Let us first consider the corrections arising from two insertions of Vs (at second order time- 
independent perturbation theory). As the correction arises from potential momenta the soft 
renormalization scale has to be used. We find that the relative corrections are at most 2.8% 
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FIG. 9. Variation of the NNLL cross section for a) the value of the strong coupling, b) the top 
quark width, and c) the inclusion of a Standard Model (SM) Higgs boson. Changes relative to the 
central value (solid lines) are shown by dashed red lines. In c) there are two solid lines, the lower 
black line is the decoupling limit for the Higgs boson, and the upper blue line is for a SM Higgs 
with mass mu = 115 GeV. 



(for v = 0.1 and close to the peak). In Ref. [Q the corrections to the square of the .S-wave 
quarkonium wavefunction at the origin arising from the emission and reabsorption of an 
ultrasoft gluon were determined. Naturally the ultrasoft scale has to be used in this case. 
These corrections can be taken as an estimate for the ultrasoft corrections to the normal- 
ization of the total cross section. For the ground state (n = 1) one finds that this ultrasoft 
correction amounts to about 2% for v ss a,. This is also consistent with our error estimate. 



VIII. PHENOMENOLOGICAL CONSEQUENCES 

Our result for the total cross section confirms the previous result in the fixed order ap- 
proach that a top quark mass measurement with theoretical uncertainties at the level of 
100 — 200 MeV or better is feasible [|17|- In principle, the prospects are now even better, 
because the small normalization uncertainties of the renormalization group improved cal- 
culation have eliminated the danger that normalization uncertainties could feed into the 
top mass measurement through beam effects caused by initial state radiation, beam en- 
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ergy smearing, and beam-strahlung. On the other hand, it is also natural to reconsider the 
threshold scan as an instrument to carry out precise measurements of other top quark pa- 
rameters such as the strong coupling a s , the top Yukawa coupling y t or the total top width 
Y t - ideas that have not been seriously considered since the large NNLO corrections to the 
normalization of the cross section were discovered. In fact, in view of our renormalization 
group improved NNLL calculation, older studies assessing how precise these parameters can 
be measured (see e.g. Refs. |HH]), appear more realistic than conclusions based on the NNLO 
predictions. 

Here we discuss the sensitivity of the total cross section a^ t to changes in the input 
parameters a™ l=5 (m z ), Y t , and y t (for a light Higgs). Our discussion is restricted to the 
Standard Model. However, we note that low energy supersymmetry does not lead to con- 
siderable changes in the QCD dynamics p9" |, even though it can lead to modifications of 



the value of electroweak parameters that can be measured at the top threshold such as T t 
and y t . For the implementation of the Higgs effects we use the prescription of Ref. [S0|. For 
Higgs mass m H , the Higgs potential 



y 2 t i 

is accounted for exactly through Eq. fl3"5|), and the one-loop short- distance corrections to the 



Vm(p,q) = -f ry- 2 (52) 



current O p \ that depend on y t are included in the matching condition Ci(l) [|60] 



+ ;(ir7j7(l- 5 + rf) cos -.(^ 



(53) 



where r = mnjvn t . Assuming that m# has been measured elsewhere, the question to address 
from the top threshold cross-section is how well y t can be measured. 

In Figs. |9] we have displayed the variation of the NNLL cross section cr7dt as a function 
of y/s for different choices of the input parameters a s (mz), Y t and y t . In all figures we have 
chosen M t is = 175 GeV, T t = 1.43 GeV, a s (m z ) = 0.118, y t = and v = 0.15, unless stated 
otherwise. The upper left panel shows the cross section for 0.116 (lower dashed red line), 
0.118 (solid black line), and 0.120 (upper dashed red line). At the peak one finds a ±2.7% 
variation when varying a s (mz) by ±0.002. The upper right panel shows the cross section 
for Y t = 1.43 GeV (solid black line). The dashed red lines correspond to variations of Y t 
by ±10%, where for a smaller width the peak becomes more pronounced. At the peak one 
finds a (—2.3%, ±2.6%) variation of the cross section when varying Y t by ±5%. Finally, the 
lower panel shows the cross section for zero Yukawa coupling (solid black line), the Standard 
Model value for the Yukawa coupling (solid blue line), and a ±20% variation of the coupling 
with respect to the SM value (upper/lower dashed red lines). The Higgs mass was chosen to 
be m# = 115 GeV. At the peak one finds a (±3.3%, —2.6%) variation when varying yf M by 
±20%. For ran = 130 and 150 GeV the corresponding variations are (±2.9%, —2.3%) and 
(±2.5%, -2.0%) respectively. 

Thus, given the remaining 3% uncertainty from higher order QCD effects, measurements 
with uncertainties of 5a s {m z ) ~ 0.002, 5Y t /Y t ~ 5% and Sy t /y t ~ 20% appear feasible. More 
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accurate extractions might be possible if better use is made of the smaller scale dependence 
above and below the IS peak. A detailed simulation study based on the NNLL results will 
be necessary to more realistically assess the uncertainties in measurements of a s , T t and y t . 
Such a simulation is beyond the scope of this work. 

Finally, we comment on other prospects for measurements of a s , T t , and yt- Extractions 
of a s from many processes currently have a combined ±0.002-0.003 uncertainty |6T1 . |62| , but 



in the very long term a 1% uncertainty may be achievable J55, 51|. Measurements of the 



forward-backward asymmetry and momentum distribution for top quarks near threshold 
may help to reduce the uncertainty in an extraction of T t [R EJ . The top Yukawa coupling 



can also be measured above threshold from e + e~ — > tih p5| . However, at ^fs = 500 GeV 
the phase space for this process is restricted, so precise extractions above threshold require 
a 0.8-1 TeV linear collider and high luminosity. 

IX. SUMMARY AND OUTLOOK 

In this work we have studied the impact of the summation of QCD logarithms of ratios of 

the scales m t , m t v and m t v 2 on the total top pair production cross section in the kinematic 

regime close to threshold. Using the recently developed effective field theory vNRQCD 

22| p5[ the logarithms of velocity are summed in the Wilson coefficients of the vNRQCD 



operators based on renormalization group equations that are determined from the anomalous 
dimensions of the nonrelativistic operators. All NNLL terms were accounted for, except for 
the Wilson coefficient c\ of the dimension three production current for which only NLL results 
are completely available [p5" |. Based on partial contributions to the three loop anomalous 



dimension that are currently known, the modification of c\ at NNLL order was estimated to 
be at the 1% level. We find that the size of the NNLL corrections and the variation of the 
NNLL cross section for different choices of the scaling parameter v is an order of magnitude 



smaller than the results of earlier NNLO calculations ||17|| , where the logarithmic terms were 
not summed. We conclude that a conservative estimate of the remaining theoretical QCD 
uncertainty of the normalization of the cross section is ±3%. The already excellent prospects 



of a top quark mass measurement at the level of 100-200 MeV or better [0 , are in principle 



further improved because the danger that normalization uncertainties could feed into the 
top mass measurement through beam smearing effects is eliminated. Also, measurements of 
a s , T t and y t from a threshold scan with theoretical uncertainties of about 5a s (mz) ~ 0.002, 
ST t /T t ~ 5% and 6y t /y t ~ 20% appear feasible at this stage. However, realistic simulation 
studies are still necessary to obtain more definite numbers. 

At the level of uncertainty of our computation of the total cross section there are addi- 
tional subleading effects which should be reconsidered. Many small effects were previously 
considered irrelevant for the total cross section due to the large size of the QCD corrections 
found at NNLO. Here, let us just mention electroweak effects. The dominant contribution 
was included and comes from the top width, which is the reason that the threshold top 
cross section is a smooth lineshape and can be calculated for all center of mass energies. 
At subleading order a number of additional effects need to be taken into account, such as 
the so-called non-factorizable corrections |3TI], electroweak box and triangle diagrams K> 



W- width effects [0 , and single- or non-resonant background processes that lead to the same 



final state as top pair production and phase space corrections |p.5| . These corrections lead 
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to effects at the percent level and are expected to be of the same size as the remaining QCD 
uncertainties. 
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APPENDIX A: CONSTANTS AND RUNNING COUPLINGS 



1. Formulae for constants 



In this appendix expressions for the constants that appear in the effective Coulomb 
potential V c (p,q) are collected and the two loop matching for Ci(l) is discussed. The (3- 
functions and constant coefficients which appear in V c (p, q) in Eq. (|TT|) include 
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where the a^ term was computed in Ref . [0 , and ni = 5 is the number of light fermions. 

For the matching onto the coefficient ci(l) of the non-relativistic current O p i in Eq. (|T1]). 
the one-loop contribution is well known. The two-loop computation involves computing the 



difference of graphs in full QCD and the effective theory as shown in Fig. 10. We find 



d(l) = l 



2CV 



TX 



a s (m) + a 2 s {m) 



C F ( 



In 2 



12 



25 
24 



7T- 



+ C A C F (ln2-l + 



K 



(A2) 



K — Lj p 

+C F Tp 



1/39 A x 4, 35 



4/11 

9^ _1 



CpTpm 



— c A c F 

11 



1 /151 13 n 8, 



179 
~T2 



9tt 2 



At order a s (m) 2 this matching result is scheme dependent and depends both on the sub- 
traction scheme as well as on the definition of operators in the effective theory. For this 
reason the result in Eq. ( |A"2]) differs from the earliest results in Refs. [p8| ,f2~9H for Ci(l). In 
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FIG. 10. Difference of full QCD |11] and order a 2 v° EFT graphs which gives the two loop 
matching for ci(l). The x denotes an insertion of the p 4 /(8m 3 ) operator, and graphs with this 
insertion on a different propagator are understood. 



Refs. p8| , ^9| the potential operators in Eq. fll2|) were denned with dependence on d = 4 — 2e 
in such a way that only the hard part of the two loop graphs in the threshold expansion 
contributed. The result in Eq. ( \K2ft uses the MS scheme with all effective theory operators 
defined on-shell in d = 4 space-time dimensions. The full theory result was extracted from 



Ref. ||28|| ) while the EFT graphs were computed explicitly. We have obtained the same result 
by using the "direct matching" of Ref. |]£)8|| , where one matches results for cross sections in 
the full and effective theories rather than matching the Green functions. 



2. Formulae for running couplings 



In this section the formulae for the Wilson coefficients of operators needed for our NNLL 
cross section are given. The Wilson coefficients for the potential operators in Eqs. ([TT|) and 
(|j~2|) were computed in Refs. [|23],[25|.f26|, and due to mixing depend on the one- loop running 
of the HQET 1/m 2 terms computed in Ref. |6JJ. For completeness we summarize the results 
for the color singlet coefficients: 
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FIG. 11. Examples of the type of graphs that contribute to the three- loop anomalous dimension 
for c\. 5V2r are one loop counterterms. 



Vi s) ( 



v 



-A7iCFza s (m) 



where 



32ttCfCa 



aAmv) 



64tcC f C a 

®s{m) H — \n(w) a s {m) 

6p 



a s (m) 



w 



a s {mv 2 ) 
a s {mv) 



(A3) 



(A4) 



Next consider the Wilson coefficients of the currents in Eqs. (14) and (|15|) . The LL 
solutions are C\{v) = c^v) = 1, while the solution for C2{y) is given in Eq. fllpp. The NLL 
solution for C\{y) is obtained by integrating the two-loop anomalous dimension [ [22] 
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with one-loop matching. The analytic solution of Eq. ( |A5| ) can be found in Ref. |25| . 

At NNLL we use the two-loop matching in Eq. ( |A2| ) , however the complete three- loop 
anomalous dimension is currently unknown. Part of the unknown contributions are due to 
the mixing in Eq. (|A5|) . Solving for the full NNLL c\(v) requires the NLL values of V c (z/), 
V2(^), V r (z/) and V s (^), and the NNLL value of Vfc(z/), which in turn require knowing the 
running of 1/m 2 operators in the HQET Lagrangian at two-loops. Note that at one-loop, the 
calculation of the full mixing matrix in Ref. |BI| already involved 75 diagrams. In addition 
there are direct contributions to the three loop anomalous dimension from three-loop graphs 
with divergences which require a C\ counterterm. By power counting we find that there are 
no contributions from purely potential three loop diagrams. However, diagrams with soft 
and ultrasoft gluons such as those shown in Fig. [11] may contribute. Graphs with ultrasoft 
gluons are expected to generate an anomalous dimension of the form 
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where A is (currently) an unknown number. A complete calculation of A may not be that 
difficult, and is of interest since the graphs with soft gluons give contributions which are 
parametrically smaller by a factor a s (mv) / a s (mv 2 ) . 

Interestingly, we note that the cancellation of sub divergences ln(p)/e and \n(E)/e between 
Fig. |TT| b and the counterterm graph in Fig. |TT| c requires the velocity renormalization group 
relation jjl v = [i 2 s /m for consistency. In Coulomb gauge 
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It is easy to see that the sub divergences of the form ln(- • -)/e cancel completely only after 
using the equations of motion E = p 2 /m and pu = p"s/m. The difference of #/e terms in 
Eq. ( |A7| ) will contribute to Eq. (|A6[). If only a single /i were used then one would be left 
with a ln(/z/m) term contributing to the anomalous dimension which is obviously wrong. 

To numerically estimate the effect of the three loop anomalous dimension on Ci(v) we 
consider several types of contributions: 

1. At NLL the running of V^ s '(v) is determined by the two-loop (3- function, so a subset 
of terms in the anomalous dimension are determined by increasing the order of the 
running used for a s . Using one-loop running for z and w in Eq. 



gives cf NLL (z/ 



0.15) = 0.8970, whereas two-loop running gives c\ 



NNLL/ 



0.15) = 0.8941. 



(»), 



2. Another estimate comes from partial knowledge of the NLL results for V 2s (^) and 



the NNLL result for V, 



W, 



v 



which enter in Eq. 



The matching for these poten- 



tials is known, however the complete anomalous dimensions are not. The one- loop 
matching values are |2|, [H| V : 



«/ 



= 241/135 a 2 s (m) and V s (s) (l) = 112/27 a 2 s (m). 
The recently computed two-loop matching [71| for the l/|k| potential gives Vfc(l) = 
— (289 + 4896 ln2)/(3247r) a^(m). Furthermore, the QED anomalous dimensions for 



these potentials are known p3|, and directly give a majority of the purely potential 
contributions to the QCD anomalous dimensions (in the notation of Ref. |33[] we in- 
clude the terms p ccc , p cc2 , p c22 , p c k, and p^)- Using these results as approximations for 
v) and V;[ (v), and then including these contributions in Eq. (|A5|) gives 



yWNLL 
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0.15) 
0.15) 



0.8991. The result without the two-loop matching from Ref. [[7Bj is 
0.8915. 



3. Finally, adding a contribution of the form in Eq. 



possibility of a large coefficient gives the range q" 
(using three loop running for a s ). 



NNLL/ 



with A = ±5 to account for the 
= 0.15) = 0.8903 to 0.8969 



Adding the magnitude of these three anomalous dimension contributions linearly, we find 
a 1.2% change in c\{y = 0.15) relative to the NNLL value which includes only the two- 
loop matching for c\. We therefore do not expect that the missing three loop anomalous 
dimension for c\ will shift the value of R v by more than an amount at the 2% level. 
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